Disease spreading in populations of moving agents 
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We study the effect of motion on disease spreading in a system of random walkers wliich addition- 
ally perform long-distance jumps. A small percentage of jumps in the agent motion is sufficient to 
destroy the local correlations and to produce a large drop in the epidemic threshold, well explained 
in terms of a mean-field approximation. This effect is similar to the crossover found in static small- 
world networks, and can be furthermore linked to the structural properties of the dynamical network 
of agent interactions. 

PACS numbers: 89.75.Hc, 89.75.-k,87.23.Ge 
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Many information/communication and social systems 
can be modeled as complex networks [H, 0, Q • One of the 
main reasons for studying such networks is to understand 
the mechanisms by which information, rumors and dis- 
eases spread over them. Recent works have pointed out 
the importance of incorporating the peculiar topology of 
the underlying network in the theoretical description of 
disease spreading 0j S ■ Epidemic models are in fact 
heavily affected by the connectivity patterns characteriz- 
ing the population in which the infective agent spreads. 
Both the nature of the final state, and the dynamics 
of the disease process, strongly depend on the coupling 
topology. Specifically, spreading occurs faster in small- 
world systems, i.e. in networks with shorter characteristic 
path lenghts [tI, fij. Moreover, the epidemic threshold is 
affected by the properties of the degree distribution P{k). 
For instance, the divergence of the second-order moment 
of P{k) leads, in uncorrelated scale-free networks, to the 
surprising result of the absence of an epidemic thresh- 
old and its associated critical behavior [l^, [llj . This 
implies that scale-free networks are prone to the spread- 
ing of infections at whatever spreading rate the epidemic 
agents possess. 

Most of the results present in the literature so far refer 
to cases where the disease spreading takes place over a 
wiring topology that is static, i.e. the underlying network 
is fixed in time, or grown, once forever. A more realistic 
possibility is to consider the networks themselves as dy- 
namical entities. This means that the topology is allowed 
to evolve and adapt in time, driven by some external fac- 
tors or by the very same spreading process. For instance, 
Refs.[T2l, [Tsl have considered disease spreading on adap- 
tive networks in which the susceptible agents have per- 
ception of the risk of infection, and are able to avoid 
contact with infected agents by rewiring their network 
connections. In this Letter we study disease spreading on 
a system of mobile agents. The agents are random walk- 
ers which can additionally perform long-distance jumps, 
and are only able to interact with agents falling within 
a given interaction radius apart from them. Hence, the 
interaction network between individuals is a dynamical 



one, because the links evolve in time according to the 
agent movement. The focus of our work is on the influ- 
ence of the motion on the disease spreading. With the 
aim of large scale simulations, we will show that the mo- 
tion, usually neglected in epidemic models, has instead a 
profound effect on the dynamics of the spreading and on 
the epidemic threshold. In particular, a small number of 
long-distance jumps in the agent motion is sufficient to 
produce a large drop in the epidemic threshold, as that 
observed in static small- world networks [8,]. The case of 
infecting moving individuals has been considered only in 
few other works [H, [H, [H, [O, [H, [H , although the ef- 
fects on the epidemic threshold have never been studied 
in detail before. This is a quite important issue both 
in social and in artificial networks. E.g., recently the 
analogy with epidemic spreading has been exploited to 
propose routin g al gorithms in highly mobile networks of 
computers [20II2H 

We consider a system of N identical agents indepen- 
dently moving in a two-dimensional cell of linear size 
D, with periodic boundary conditions. Fixing the value 
of D is equivalent to fix the agent density p = N/D^. 
The agents are represented as point particles, and their 
positions and velocities at time t are indicated as ri{t) 
and Vi(t) = {vi{t) cos di{t),Vi{t) sin 6i{t)), i — 1,...,N. 
We further impose that the agents move with a veloc- 
ity modulus which is constant in time and equal for all 
the agents, i.e. Vi{t) — v,Wi — l,..,iV and \/t. At time 
t = the N particles were distributed at random. At 
each time step, the agents change stochastically the di- 
rection angles 9i (t) . The positions and the orientations of 
the particles are thus updated according to the following 
rule: 



0^{t) = 



(1) 



where are N independent identically distributed ran- 
dom variables chosen at each time with uniform probabil- 
ity in the interval [— tt, tt]. In addition, to include the pos- 
sibility that agents can move through the bidimensional 
world with time scales much shorter than those related 
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to disease, as in the case of infected individuals travel- 
ling by flights [13], we consider that agents can perform 
long-distance jumps. This is accounted for by defining 
a parameter, pj S [0,1], that quantifies the probability 
for an agent to perform a jump into a completely ran- 
dom position. In summary, at each time step, each agent 
evolves: following Eqs. ([1]), with a probability 1 — pj, 
or performing a jump, with probability pj. In the lat- 
ter case the position of the agent is updated into a new 
position chosen at random in the cell. Models with dif- 
ferent jumping rules [l^, non-indentical and interacting 
agents, have also been considered, and the results will be 
reported elsewhere. Finally, the main parameters con- 
trolling the moving agents in our model are: p, v and 

Pj- 

Among the possible mechanisms of disease spreadings 
d, H, 0] , we focus on the SIR model, that divides the TV 
agents into three disjoint groups: susceptible (S), infec- 
tive (I) and recovered (R). We indicate as Ns{t), Nj{t) 
and Nii(t), respectively, the number of agents in the 
three groups at time t, with the total number of agents 
Ns{t) + Ni{t) + N]i{t) = N being constant in time. A 
small number of agents is set in the infective state at 
t = as the seed of the infection, while all the others 
start from the susceptible state. The process through 
which the disease spreads can be summarized as follow- 
ing. An interaction radius r is fixed (r = 1 in all our 
calculations), and this defines the interaction network: 
at each time step t each agent interacts only with those 
agents located within a neighborhood of radius r. For a 
given susceptible agent, the probability of being infected 
increases with the number of infected individuals in the 
neighborhood. More precisely, if an agent is in the S 
state at time t, and exactly one of its neighbors is in the 
/ state, then it moves into the / state with probability A 
and stays in the S state with probability 1 — A. If Nj^ 
is the number of infected individuals in the neighbor- 
hood of the agent, then its probability of being infected 
is 1 — (1 — A)^^"- . In addition to this, each infected agent 
can move into the R state with probability /i, and then 
cannot catch the disease anymore. This sets the average 
duration time of the infection: t = -. 

In our model we implement the motion rules and we 
update, at each time step, the disease state of every 
agent. The model is simulated for a number of time 
steps sufficiently high to ensure that, at the end, there 
are no more infected individuals in the population. Dur- 
ing a simulation, the number Nj{t) of infected grows up, 
reaches a peak value, and then decreases. Typical cases 
are shown in Fig. [H where it can be noticed that a larger 
value of pj increases the spread of the infection. In fact, 
the peak of infected individuals is higher in the pres- 
ence of a larger probability of jumping. We have verified 
that this is also true for the total number of individu- 
als which have contracted the disease at the end of the 
process. Both the two issues have important practical 
consequences, since on one hand the disease involves a 
higher percentage of the population, and, on the other 
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FIG. 1: (Color online). Number of infected individuals as 
a function of time t. We have considered a system with 
N = 1000, p = 1, V ^ 0.1, A = 0.1, n = 0.05 and differ- 
ent values of Pj. The system is started with 1% of the agents 
set in the infective states. Results are averages over 100 runs. 
The dashed line is the result in the homogeneous mixing ap- 
proximation of Eqs. ((2J 

hand, it requires more resources to deal with a higher 
peak of infected individuals. 

The behavior for large pj and/or large v can be inter- 
preted in terms of a mean-field approximation. In fact, in 
such limit, we expect that the spatial correlations in the 
disease states are destroyed by the agent motion, and we 
assume that the homogeneous mixing (HM) hypothesis 
is valid. Under this hypothesis, all the individuals have 
the same probability of contacting any other individual 
(i.e., the population mixes at random) and 
the equations for the system of infectious mobile agents 
read: 

i{t+\) = i{t) + s{t) [1 - {1 - Xi{t)f] - fii{t) 

r{t + l) = r{t)+Xt(t) (2) 

s{t+l) ^ p-i{t+l)-r{t + l) 

where s{t) = Ns{t)/D^, i(t) = Ni{t)/D^, r(t) = 
Nji{t)/D^, are respectively the densities of susceptible, 
infected and recovered individuals at time i, and a = ■nr^ . 
The third equation is simply derived from the conserva- 
tion of the number of agents. The second equation indi- 
cates that the increase of recovered individuals at time 
t-f 1 is proportional to the number of infected individuals 
which get recovered, i.e. to Xi(t). The first equation can 
be derived by taking into account that the density of in- 
fected individuals at time t -I- 1 is decreased by Xi(t) and 
increased by the density of susceptibles catching the dis- 
ease. This last term is proportional to s{t) times a conta- 
gion probability pcont ■ The contagion probability is given 
by Pcont = 1 - Pcont, where Pcont represents the proba- 
bility of not being infected, pcont is the probability that 
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an agent will not being infected by none of its neighbors, 
i-e. pcont ~ ~ -^nni^) ' with a representing the area 



in which each agent may sense other infective individuals, 
that in our case is equal to 7rr^ . The number of infective 
agents as a function of the time in the HM approximation 
computed from Eqs. with A = 0.1 and ^ — 0.05, is re- 
ported as a dashed line in Fig.[T] As expected, the curves 
for the model approach the dashed line when pj 1. 
The mean field approach gives us also information on the 
epidemic threshold. In fact, for small i{t) we can approx- 
imate (1 — Ai(t))° ~ 1 — aXi{t), and we get from Eqs. ^ 
the iterative rule i{t + 1) = i{t) + irr^ \s{t)i{t) — iii{t). 
This allows to distinguish two cases. In fact, by assum- 
ing s(0) ~ p, we get > z(0) when a 



while < i{0) when a = A/^ < 
HM hypothesis we derive a critical threshold: 



7Tr-p 



Thus, under the 
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(3) 



Hence, when a < ac, the number of infected decrease 
monotonically, while for cr > an epidemic outbreak 
occurs. Also notice that CTc = as found in Erdos 

and Renyi random graphs [1, U where < fc > is the 
average number of first neighbours of an agent, that in 
our case is equal to pnr'^. 




FIG. 2: (Color online). Epidemic threshold as a function of 
the system density p. We have considered A'' = 1000, v = 0.1, 
and p — 0.05. Different curves correspond to different values 
of Pj . Results are averages over 100 runs. The dashed line is 
the prediction in the homogeneous mixing approximation. 

In Fig. [5] we report the epidemic thresholds ctc com- 
puted numerically for the model of infective agents mov- 
ing with different values of pj . For a given density p, We 
observe that Cc is a decreasing function of the density p. 
Moreover, for a given value of p, the threshold decreases 
with the jumping probability pj . In the same figure we re- 
port for comparison the prediction of Eq. ([31), as a dashed 



line. We notice that the homogeneous mixing approxi- 
mation becomes more and more accurate when pj tends 
to 1. The convergence to the HM threshold of Eq. ^ as 
a function of pj is rather fast. For instance, already at 
Pj =0.1, the threshold in the model is, for any value of p 
reported, practically indistinguishable from the HM one. 
In Fig. [2] we have considered a fixed velocity v = 0.1. 
We have also studied <Tc as a function of w, at a fixed 
density and for different values of pj . We observed that 
for small values of v, the epidemic threshold tends to the 
prediction of the HM when pj tends to 1, while, for large 
enough values of v, the epidemic threshold is consistent 
with the homogenous mixing one, independently of pj. 
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FIG. 3: (Color online). Scaled epidemic threshold 
'^<:{pj)/o'c{Pj = 0) as a function of pj. We have considered 
V — 0.1 and p = 0.05. Different curves correspond to dif- 
ferent values of the density p. Results are averages over 100 
runs. 



Finally we investigate in more details the effects of pj 
on the epidemic threshold. In Fig.[3]we report, as a func- 
tion of Pj , the value of the threshold CTc {pj ) , normalized 
by ac{pj — 0), for different values of the density p. We 
observe a rapid drop in the curves (note the logarithmic 
scale for pj ) , meaning that a small number of long-range 
jumps produce a large decrease in the epidemic thresh- 
old. The plateau observed for pj larger than 10^^ implies 
that in order to have a significant change in the epidemic 
threshold (and in the disease incidence) in our model 
of moving agents, the jumping probability has to be ex- 
tremely small. This is in line with other results, stressing 
the role of the large scale properties of the airline trans- 
portation networks in determining the global diffusion 
pattern of emerging diseases \v?\ , and can have important 
implications in the immunization of real communication 
networks [23 | . Different curves in the figure correspond to 
different densities. As p decreases, the drop in the curve 
occurs for smaller and smaller values of pj, suggesting 
that no finite critical value of pj can be determined this 
way. This behavior is similar to the crossover observed in 
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the characteristic path length of small-world networks as 
a function of the rewiring probability 0, HI] . We notice, 
however, that in our case the effect is due to the agent 
movement and not to the rewiring of static links. 
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FIG. 4: Characteristic path length Lg (circles) and clustering 
coefficient Cg (squares) as a function of Pj (A'^ = 1000, p = 1, 
u = 0.1). Results are averages over 10 realizations. 

The observed behaviour can be related to the topo- 



logical properties of the underlying dynamical network. 
For such purpose we define an effective adjacency matrix 
Gr{t) — {gij{t)} taking into account that each infected 
individual may infect other individuals during the aver- 
age duration of the infection, i.e. during t — — simula- 
tion steps Let A{t) be the adjacency matrix at time t 
defined so that [t) = 1 if the j*'* agent is within the in- 
teraction radius of the z*'* agent at time t, and {t) — 
otherwise. We set gij{t) ~ 1, if at least for one t' , with 
t' = t,t — 1, ■ ■ ■ ,t — T + 1, it is verified that aij{t') — 1. 
Otherwise we set gij{t) — 0. In Fig.|4]we report the char- 
acteristic path length, Lg, and the clustering coefficient, 
Cg of matrix G{t) as a function of pj. The behavior we 
observe is in all similar to that found in small-world net- 
works for increasing rewiring probability Q . Notice that 
the drop in Lq occurs at the same value of pj at which 
we have found the drop in the epidemic threshold. 



In summary, we have considered a system of infectious 
mobile agents to study the effects of long-range moves on 
the disease spreading. Our results indicate that the in- 
terplay between dynamics and topology can have impor- 
tant consequences for the global spreading of infectious 
diseases in systems of mobile agents, and in the related 
applications such as the forecast of epidemic spreading 
and the development of wireless routing strategies. 
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